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Abstract: We simulate five different systems belonging to the universality class of the 
gauge dual of three-dimensional random percolation to study the underlying effective string 
theory at finite temperature. All the data for the finite temperature string tension, when 
expressed by means of adimensional variables, are nicely described by a unique scaling 
function. We calculate the first few terms of the string tension up to order T 6 and compare 
to different theoretical predictions. We obtain unambiguous evidence that the coefficients 
of T 2 and T terms coincide with those of the Nambu-Goto string, as expected, while the 
T 6 term strongly differs and is characteristic of the universality class of this specific gauge 
theory. 
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1. Introduction 

The possibility of describing the long-distance dynamics of strong interactions in the con- 
fining phase by an effective string theory is a fascinating, many years old conjecture which 
dates from before the formulation of QCD In gauge theories it is based on the very in- 
tuitive assumption that the colour flux connecting a pair of distant quarks is concentrated, 
in the confining phase, inside a thin flux tube, which then generates the linear rising of 
the confining potential. According to the common lore, this thin flux tube should behave, 
when the quarks are pulled very far apart, as a free vibrating string ||]. 

The string-like nature of the flux tube is particularly evident in the strong coupling 
region, where the vacuum expectation value of large Wilson loops is given by a sum over 
certain lattice surfaces which can be considered as the world-sheets of the underlying con- 
fining string. At the roughening point || this sum diverges and the colour flux tube of 
whatever three-dimensional or four-dimensional lattice gauge theory undergoes a transition 
towards a rough phase. It is widely believed that such a phase transition of the flux tube 
belongs to the Kosterlitz-Thouless universality class Accordingly, the renormalisation 
group equations imply that the effective string action S describing the dynamics of the flux 
tube in the whole rough phase (the one connected with the continuum limit) flows at large 
scales towards a massless free field theory. Thus, for large enough inter-quark separations 
it is not necessary to know explicitly the specific form of the effective string action S, but 
only its infrared limit 

S[h] = S cl + S [h] + . . . , (1.1) 

where the classical action S c i describes the usual perimeter-area term, h denotes the two- 
dimensional bosonic fields /ij(£i, £2) with i = 1,2, . . . ,d— 2 which describe the transverse 



- 1 - 



displacements of the string with respect the configuration of minimal energy, £1 , £2 are the 
coordinates on the world-sheet and Sq[Ji] is the Gaussian action 

So[h] = \Jd 2 id a h l (i 1 ^ 2 )d a h i (i 1 ^ 2 ) (a = 1,2; i = 1, 2, . . . ,d - 2) . (1.2) 

In this IR approximation the effective string is known as the free bosonic string. The 
ensuing universal string fluctuation effects || were first unambiguously observed many 
years ago in the 7Li gauge theory in three dimensions |7|, [8| . 

In order to study the first perturbative corrections to the IR limit it has been often as- 
sumed, for the sake of simplicity, that the effective string action is the Nambu-Goto action, 
i.e. the one proportional to the world-sheet area. Expanding in the natural dimensionless 
parameter l/(aA), where a is the string tension and A the area of the minimal surface 
bounded by the Wilson loop, one can write 

S[h] = S d +S [h} + ^ J d 2 ^ [[djiid^f - 2d a h i dPh i d a h j dphi]+0 (7^2) ■ (°) 

Also in this case the first numerical analysis has been performed in a 3D Z2 gauge model 
|| . More recently, high precision numerical simulations in SU(N) gauge theories confirmed 
these effects in the static quark potential flC| , |ll]]. Mismatches between the observed 
spectrum of the low-lying string states with fixed ends and the predictions of the free 
bosonic and the Nambu-Goto strings have been repeatedly reported |l^, [L3|, [L4|, 15|. Recent 
studies demonstrated that these mismatches are gradually disappearing at larger distances 
H,|7|. 

Closed strings wrapping around a compact dimension [18| or strings at finite temper- 
ature [19] have also been considered. In this case a remarkable agreement between the 
observed data and the Nambu-Goto predictions has been reported. From a theoretical 
point of view the reasons of this agreement can be understood, at least in part, resorting to 
a systematic expansion of the most general form of the effective string action S[h] in terms 
of h a {^) and its derivatives [JO], [23]. The outcome of these studies can be conveniently sum- 
marised in some general properties of the first few terms of the low temperature expansion 
of the string tension 

a(T) =a -(d- 2)^T 2 + £ s n T n . (1.4) 

n>3 

The second term on the right hand side is the low temperature analogue of the Liischer term 
of the inter-quark potential |pij| . It is a characteristic quantum effect of the IR free string 
limit (hi) and it is expected to be independent of the interaction terms of the effective 
theory. On the side of the gauge theory it is more than universal, in the sense that it does 
not depend on the nature of the gauge group. 

As a consequence of a certain open-closed string duality 1 2C ] it was shown that for any 
number of space-time dimensions S3 = and that in three dimensions S4 is again a more 
than universal coefficient which can be evaluated in various ways [22, ^] and coincides 
with the Nambu-Goto value s± G 

SA = sf G = -(d-2) 2 -^. (1.5) 
72 cr 
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A different approach to effective string theory [gj ] leads to similar conclusions [25, ^] (i.e. 
S3 = 0, S2 and S4 more than universal), but for all values of d. 

In spite of the remarkable agreement of the first few terms of the Nambu-Goto expan- 
sion with the numerical results, theoretical reasons indicate that the Nambu-Goto string 
is a sick theory and cannot describe the effective confining string to all orders in T: de- 
pending on the quantisation method, one finds either the breaking of rotational invariance 
or appearance of the conformal Liouville mode, in contrast with the assumption that the 
only physical degrees of freedom of the effective string are the transverse modes. 

Numerical experiments lead to similar conclusions, showing that different gauge theo- 



ries are described, at least at short distance, by different effective strings |14|, pL5 1 , even if, 
so far, the order of the first term deviating from the more than universal behaviour in the 
power expansion of cr(T) has not been determined. In this paper we find the order of such 
a term by evaluating in a particularly simple model, the gauge dual of random percolation 
in three dimensions, the coefficients s n up to n = 6 order. We find that sq strongly deviates 
from the value predicted by the Nambu-Goto model. 

We performed five different kinds of high-precision numerical experiments by varying 
the implementation of the percolation, the lattice spacing, the temporal extent of the lattice 
and the type of lattice. All these variations should keep the system in the same universality 
class. Indeed, as expected, all the collected data agree with the more than universal values 
of S2 and S4 and lead to 

(See Table § for more details). The vanishing of S5 suggests that the high temperature 
expansion is even in T, like in Nambu-Goto case. Notice however that the value we find 
for sq for the gauge dual of percolation is very different from the corresponding coefficient 
of Nambu-Goto string, which is negative: Sq G = — (d — 2) 3 ^a' 2 ' P renmmarv results have 



been presented in [27, 28]. 
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2. Polyakov loops 

We focused on the behaviour of the Polyakov-Polyakov correlation function at finite tem- 
perature in a (2+l)-dimensional system. The lattice is a L 2 x i slice with periodic boundary 
conditions, with L large enough to represent the spatial extent and £ = the inverse tem- 
perature. We considered a pair of Polyakov loops orthogonal to the spatial direction and 
at a distance of r lattice spacings a; the (connected) correlation function in this case is 
denoted by (P(0)P*(r) ). 



At the free string or leading order (LO) approximation (1.1) the functional form of 



this correlator in the effective string picture was calculated in different contexts. In lattice 



gauge theory it was first derived in [29|, leading to 

-air— ill 



(P(0)P*(r)) LO K e ^ T)d _ 2 , (2.1) 
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where the Dedekind rj function is defined as 



1 , if 
V{r) = q^H(l-^), r=-, q ^ e 2 ™ . (2.2) 

n=l 

Within this approximation the temperature-dependent string tension, defined as the coef- 
ficient of the linear part of the confining potential, i. e. 

a(T) = - lim -1 log<P(0)P»> , (2.3) 

i — >oo rl 

turns out to be 

a (T) = a -( d -2)^T 2 = a -(d- 2) -T 2 + 0(T 4 ) , (2.4) 
b o 

as expected from ( |1.4| ); cjq is the zero-temperature string tension and a = o"o + 0(T 4 ). 

At the next to the leading order (NLO) the functional form of the correlator has been 
calculated in [031 



Pimpl( e -^-^ / ( ( j-2)7r 2 £[2ii; 4 (r) + ((i-4) j E 2 2 (r)] 1 , 

(P(O)P (r))jv LO = „,_^_ 2 ( 1 + +°(l5))' ( 2 - 5 ) 



v( T ) d 2 V 1152or 3 r 

where the functions E2 and (second and fourth Eisenstein functions) are defined by: 

00 

E 2 {t) = 1-24^(7!^, (2.6) 

n=l 

00 

£!4(t) = 1 + 240 a 3 (n)q n . (2.7) 

n=l 

The functions o"j(n) here represent the sum of the i-th powers of all divisors of n. Using 
the definition (|2.3| ) one can easily verify that the parameter a is related to cr(T) and ctq 
through 

a(T) =o--T 2 - ^T 4 = a - -T 2 - -^—T 4 + 0(T 5 ) . (2.8) 
3. The model 

In this work, the model we chose as laboratory to study the effective string theory is the 
three-dimensional random percolation [30]. It can be seen as the gauge dual of a Q-state 
Potts model in the limit Q — > 1. 

It is well known that for integer Q > 1 one can formulate the gauge Potts model 
either in terms of gauge fields or in the dual version in terms of the spin variables. From 
a computational point of view the latter is much more convenient in lattice simulations. 
It is then useful to map the needed gauge invariant observables (Wilson loops or Polyakov 
correlators) into the corresponding quantities of the dual version and not to worry about 
the gauge formulation. This approach particularly well suits random percolation, as the 
direct gauge formulation is not (yet) known, but the rules to evaluate the gauge invariant 
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observables are unambiguously defined in any configuration of the system, as we shall see 
below. 

In the bond percolation model the lattice configurations are generated as follows. Each 
link of a three-dimensional lattice A is independently set to on or off according to some 
fixed probability p, which plays the role of a coupling constant. The set of on links, or 
active links, forms a graph G, whose connected components are known as clusters. When p 
exceeds a threshold value p c , depending on the nature of the lattice, an infinite, percolating 
cluster forms. 

Similarly, in the site percolation model, another possible formulation of the random 
percolation, one independently sets on or off the nodes of the lattice with a fixed probability 
p and generates a graph G by putting an active link for each pair of adjacent on nodes. 

The key ingredient to extract from the above ensemble of graphs the relevant infor- 
mation on the underlying gauge theory is the definition of the percolation counterpart of 
the Wilson operator W 7 associated with whatever closed path of the dual lattice. We set 
W 7 (G) = 1 if there is no path of G topologically linked to 7, otherwise we set W 7 (G) = 0. 
In other words, W 7 is a projector on the ensemble of graphs whose image is the subset of 
graphs not linked to 7. Therefore its vacuum expectation value (W 7 ) coincides with the 
average probability that there is no path in any cluster linked to 7. 

As in usual gauge theories, evaluating these quantities yields the main physical prop- 
erties of the model. In this way it has been shown that the percolating phase is confining. 
The string tension a and the other physical observables have the expected scaling be- 
haviour dictated by the universality class of three-dimensional percolation, therefore such 
a theory has a well-defined continuum limit Moreover it has a non-trivial glueball 

spectrum and a second-order deconfining transition at finite temperature T c with a 
ratio T c j ' yfa ~ 1.5 which turns out to be universal, i.e. it does not depend on the kind of 
lattice utilised nor on the specific percolation process considered (bond or site percolation). 

In the study of the finite size effects described by the effective string theory one could 
use in principle whatever confining gauge theory, owing to the fact that the dominant 
effects do not depend on the gauge group. The great advantage of the dual percolation 
model we study in this paper is that its simplicity allows to explore regions that are still 
inaccessible to the other gauge systems from a computational point of view. 

4. Methodology 

In this Section we describe the principal aspects of our method, based on the direct mea- 
surement of the correlator of two coplanar Polyakov loops at finite T <T C . We begin with 
the description of the lattice algorithm and proceed to discuss the kind of fits we use to 
extract the temperature-dependent string tension. 

4.1 Simulations 

We are interested in the universal properties of the effective string theory in the gauge 
dual of percolation, therefore we studied how the system responds to a variation of the 
spatial and the temporal sizes of the lattices, of the occupancy probability p, of the kind of 
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percolation (bond or site) and finally of the geometry of the lattice, considering both the 
simple cubic lattice (SC) and the body-centred cubic lattice (BCC). The set of simulations 
is listed on Table [|. The values of p are taken from and are the occupancy probabilities 



Lattice 


P 


£c = ^/aT c 


temporal sizes £ 


spatial sizes 


SC bond 


0.272380 


6 


9 -=-15 


128 


SC bond 


0.268459 


7 


10 -=-15 


128 


SC bond 


0.265615 


8 


10 -=-17 


128;194;256;320 


SC site 


0.3459514 


7 


11 -=- 17 


128 


BCC bond 


0.21113018 


3 


4 -=-10 


128 



Table 1: Relevant parameters of the simulations. 



corresponding to systems which are at the deconfining temperature aT c = l/£ c when the 
(periodic) temporal extension in units of lattice spacing a is the value t c reported in the 
Table. The simulations were made in the confined phase with a temporal extension in 
the range l c < I < 3£ c where the Polyakov-Polyakov correlator is well described by the 
NLO formula (j2.5j). The spatial size was 128 2 which was in most cases amply sufficient to 
account for the infinite volume limit. Only in two cases, namely £ = 10 and t = 11 with 
£ c = 8, we observed a non-negligible dependence on the spatial size. In those cases we 
performed further simulations on larger lattices, as indicated on the table, and extracted 
the corrected value of the string tension a using the scaling relation 

d XIL = d-cL- x ^ , (4.1) 

where = 3 is the thermal exponent of two-dimensional random percolation. In both 
cases the fit to the data was very good. 

To reach an acceptable statistics, we collected data from 10 5 configurations for each 
value of p and £. 

4.2 Algorithm 

Due to the particular nature of the random percolation model, each configuration can be 
generated independently from scratch, by simply filling an empty lattice with links (or 
sites) that are randomly switched on with a probability p. 

The tricky part is the measurement of the topological linking of the resulting graph G 
with a pair of Polyakov loops; to this end, we first choose a cylindric surface £ bounded by 
the two loops and look for the closed paths of G intersecting it and linked with one of the 
two loops. It is convenient to "clean up" the graph G — > G' , getting rid of dead ends and 
bridges between loops, as they cannot belong to the mentioned closed paths |Q. This is 
done once for the whole configuration. 

On this "minimal" configuration G' , then, the surface X is translated in all possible 
positions and the linking is measured with the technique of reconstructing each time the 
clusters in the configuration (by means of the Hoshen-Kopelmann algorithm) keeping track 
of the crossings of the loop surface, in order to detect nonzero winding numbers. 
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4.3 Fits 



The measured Polyakov-Polyakov correlators are compared with the expected behaviour 
Q2.5| ). Being this an asymptotic expression, valid in the IR limit, we fitted the data to ( |2.5| ) 
by progressively discarding the short distance correlators and taking all the values in the 
range r m i n < r < r max = 50a, with r m i n varying from the value £ indicated in the Table 
|l] to 40 lattice spacings a. The value of the fitted parameter a as a function of r m i n is 
plotted in Figure [l|. The large plateaux in the whole range of the temporal extension £ 
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Figure 1: The fitted value of a to (2.5) as a function of the minimal distance r m i n of the set of 
Polyakov-Polyakov correlators considered in the fit for the site percolation in SC lattice with I c = 7. 
The different plateaux correspond to different values of the temporal extension £. A similar plot 
for the bond percolation in the SC lattice can be found in J|7j for £ c = 6 and in J2^| for £ c = 8. 

considered show the stability of the fit which is also supported by a x 2 /dof of the order 
of 1 or less. In some cases, when £ is too close to £ c , the plateau starts at larger values 
of r m in and correspondingly the x 2 test is not good. We discarded these data from the 
further analysis. In all other cases the Polyakov-Polyakov correlator in the examined range 
of r and / is well described by the asymptotic formula ( |2.5|) . Since the latter is a result of 
the continuum, this agreement can also be interpreted as a check for the absence of finite 
lattice spacing effects at the level of our statistical accuracy. 

It is important to note that the fitted parameter a is not yet the string tension at 
zero temperature do, since ( |2.5[ ) is not an exact formula, but only takes into account the 
temperature dependence up to the order T 4 . On general grounds we expect 



a 



cto + 0(T 5 ) . (4.2) 
If it turned out that the dependence of the parameter a on T involved lower powers of T, 
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Figure 2: Plot of the fitting parameter a as a function of T 6 in numerical experiments with bond 
percolation with l c = 8. A similar plot for the case £ c = 7 can be found in pq. 



i.e. T 2 and/or T 4 , it would mean that the first two thermal corrections in fl2.8| ) were not 
universal. This question can be settled by studying the dependence on £ of the mentioned 
plateaux. In all the cases it turns out that for aT = l/l low enough the correction is 
proportional to T 6 (see for instance Figure ||). We inserted the fitted parameter a in 
J2.8|) in order to reconstruct the quantity <r(T) for the whole set of temperatures listed 
in the fourth column of Table 1. We then performed, for each line of such a Table, a 
two-parameter fit to the formula 



cr(T) = (J 



7T r, 

—T 

6 



7T- 



72a 



T 4 + 



TV 



Col 



0(T 8 



(4.3) 



The fitted parameters gq and C turn out to be stable. Their values are reported in Table 
2. Another way to analyze the data is to combine ( |4.2j ) with the observation that the term 
T 5 is absent and fit directly the parameter a to the formula a = a® + -^-^ T 6 . This way of 
analysing the data differs from the previous one for terms of the order 0(T 8 ), thus it can 
be used for a rough estimate of the systematic errors. It turns out that the evaluations 
of do coincide, within the statistical errors, with the values determined in the other way, 
while the estimates of C are about 10% larger than the values reported in Table 2. 



5. Results and conclusion 

In this paper we combined Monte Carlo simulations with different finite-size scaling tech- 
niques applied to various percolating systems. The outcome of the extensive numerical 
experiments on the gauge dual of random percolation and the analysis described in the 
previous Section is a precise determination of the string tension as a function of the tem- 
perature in a wide range of T. 
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Figure 3: Plot of the scaling variable a(T)/(Ty/ao) as a function of T/^/oq. The dashed line is 
the result of the Nambu-Goto string. 



If we plot the adimensional ratio a{T) / (T^/~oq) versus the adimensional temperature 
T I yfcjQ it turns out that all the data neatly lie on a unique universal curve as Figure 
|3| shows. This scaling behaviour indicates that the most relevant sources of systematic 
errors, including the approach to the infinite volume and the continuum limits, have been 
taken into account. The plotted quantity is expected to vanish at T c with the power law 
~ (T c — T) V2 , where v<i = | is the thermal exponent of 2D percolation. Unfortunately our 
data are not sufficiently close to T c in order to check accurately this behaviour. A similar 
scaling function has been determined for the 3D SU(2) gauge model p2|| . 



From each set of the numerical simulations described in each row of Table [T] we can 
extract three physical quantities. The first one is the coefficient C of Eq. ( |4.3| ) which 
determines the T 6 correction to the string tension. The five values of C generated by 
as many different systems (see Table |2|) remarkably coincide up to the statistical errors. 
Each set of simulations yields also a precise determination of a 2 o"o which, combined with 
the precise value of the deconfinement temperature in the same lattice units, yields the 
adimensional ratio T c /y / oo. This quantity is expected to be constant in the continuum 
limit. 

In order to extrapolate to this limit, one has to take into account the correction to 
scaling terms. The string tension in the gauge dual of percolation is expected to obey the 
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lattice i c 


— i-ialc 


L- 


a a 


X /do J 




SC bond 


6 


291(7) 


0.012612(6) 


0.15 


1.4841(4) 


SC bond 


7 


281(5) 


0.009234(5) 


1.2 


1.4866(5) 


SC bond 


8 


297(5) 


0.007059(5) 


0.4 


1.4878(5) 


SC site 


7 


307(9) 


0.009399(8) 


0.2 


1.4735(6) 


BCC bond 


3 


295(14) 


0.0474(4) 


0.8 


1.531(7) 



Table 2: The parameter C and a 2 do in the fit ( [4.3| ) and \ 2 /dof which are obtained for the 
corresponding numerical experiments listed in Table pi. The last column is the universal ratio 
T c / yf&o as obtained by combining the second and the fourth columns. 



scaling behaviour [ 30 1 



a 2 a{p) = S(p-p c ) 



2;/ 



1 



(5.1) 



i + B{p- Pc y 

where p c is the critical threshold and v and uj are the thermal and correction-to-scaling ex- 



ponents of 3D percolation (see [33] for an accurate numerical estimate of these exponents). 
Similarly, the deconfining temperature T c is expected to scale as 

1 



aT c = T(p-p c 



(5.2) 



i + c{p-p c Y 

When applied to the case of bond percolation in the SC lattice they yield S = 9.29(2) 
and T = 4.562(1), thus the extrapolated continuum limit of T c /y/aQ is estimated to be 
T/VS = 1.497(2). 

In conclusion, in this paper we extracted from various three-dimensional percolating 
systems some general information on the effective string theory describing the infrared 
properties of the confining phase of the gauge dual of percolation at finite temperature. 
We numerically evaluated the universal scaling function describing the string tension as a 
function of the temperature. We obtained clear evidence that the first two non-vanishing 
coefficients of the expansion of a in powers of T coincide with those of the Nambu-Goto 
string, while the third one strongly differs. Nonetheless this term does not depend on the 
UV cut-off nor on the specific percolation model, but is characteristic of the universality 
class of (the gauge dual of) the three-dimensional random percolation. 
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